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O I Abstract 
ON . 

' We perform a detailed numerical investigation of the approximate moduli space metric proposed 

by Diakonov and Petrov [1] for a confining model of dyons. Our findings strongly indicate that 
^ ', only for a small number of dyons at sufficiently low density this metric is positive definite - 

^ and, therefore, a valid moduli space metric ~ throughout a considerable part of configuration 

space. This poses strong limitations on results obtained by an unrestricted integration over 
collective coordinates in this model. It also indicates that strong correlations between collective 
coordinates will be essential for the physical content of a dyon model, which could be exhibited 
by a suitable simulation algorithm. 



1 Introduction 



The semiclassical approximation of path integrals around non-trivial saddle points (originating 
from the work of Callan, Dashen and Gross [21 [3] and later adapted to finite temperature by 
Gross, Pisarski and Yaffe 0]) is a prominent, but not undisputed [5l [H [7] semianalytic ap- 
proach to calculate non-perturbative effects in quantum field theories. In gauge theories at zero 
temperature the corresponding semiclassical topological objects used to be instantons [S], four- 
dimensional selfdual lumps of action density carrying one unit of topological charge. Their use 
in a semiclassical context was opened by the seminal paper by 't Hooft [9] , where the fluctuation 
determinant in the background of such objects was computed. However, the resulting distribu- 
tion with respect to the size parameter p is obviously unphysical for large instantons. Empirical 
interactions [TOl [TTl [12] have been added to the first instanton gas model [21 13] eventually leading 
to the instanton liquid model [13], where the average size is fixed to p = 1/3 fm (along with a 
density of n~^/^ = 1 fm) to roughly match phenomenological requirements. 

In pure gauge theories at finite temperature the Polyakov loop is the order parameter for con- 
finement. That there is a close relation between the Polyakov loop average and the dominating 
classical solutions has been fully realized with the discovery of calorons with non-trivial holon- 
omy by Kraan and van Baal [H] and Lee and Lu [15j . From this perspective, the previously 
known caloron solutions due to Harrington and Shepard [16j now appear as the limiting case of 
trivial holonom}0. The eigenvalues of the untraced asymptotic Polyakov loop - called holonomy 
- are responsible for the dissociation of calorons into constituents, for gauge group SU(A^). 
These are BPS monopoles with electric charges equal or opposite to their magnetic charges. We 
will refer to them as dyons (for reviews cf. [T9] ) . 

Hence, a semiclassical model of finite temperature gauge theory should be based on dyons. If the 
holonomy is related to the order parameter, all types of dyons are of equal mass in the confined 
phase but expected to split into light and heavy ones in the deconfined phase. Dyon gauge 
fields combined to intermediate size calorons do not generate confinement at trivial holonomy, 
whereas at maximally non-trivial holonomy the corresponding static potential is believed to be 
linearly rising. This has been demonstrated numerically in |20j . As size distribution in the 
confinement phase a generalization of the instanton distribution was adopted that was inspired 
by the evaluation of the fluctuation determinant in refs. |211 I22j . The suppression of heavy 
dyons due to the non-vanishing Polyakov loop above the critical temperature could be the 
mechanism behind the decrease of the topological susceptibility and the vanishing of the chiral 
condensate |23] . In contrast to this, the chiral condensate under (unphysical) periodic boundary 
conditions is non- vanishing, an effect seen in various studies |24^ [25l [26l [271 [2E] that is likely to 
reflect the presence of light dyons. 

Here we present results of a preparatory study of SU(2) gauge theory trying to obtain numerical 
insight into a dyon model recently developed by Diakonov and Petrov [Tlll91l29]. In this work, the 
role of the moduli space metric was emphasized and incorporated. The authors have generalized 
the known form of the metric of same kind dyons and of opposite kind dyons to a general metric 
valid at large distances in the parameter space of an ensemble of selfdual dyons. The problem 
of including dyons of opposite topological charge was postponed, which means that the model 
is a rather crude approximation. The authors of [1] attempted to take the effects of a mixture 
of selfdual and antiselfdual objects into account by multiplying some of their results by factors 

^For a dyon model based on these solutions see [17] , 



1 



of 2 or ^ assuming negligible interactions between the two systems of selfdual and antiselfdual 
dyons. 

The resulting moduli space metric determinant consists of Coulomb-like terms and has been 
treated analytically by rewriting the corresponding grand-canonical ensemble into an equivalent 
quantum field theory both of bosons and fermions, resembling Polyakov's famous work [30] 
showing Abelian confinement. The relations between physical quantities such as the string 
tension and the critical temperature obtained from this model are rather impressive, when 
comparing them with lattice results. However, in order to match phenomenological values, the 
model needs to be pushed to rather high densities (i.e. short distances between dyons), which 
is difficult to reconcile with the diluteness assumption of the underlying semiclassical approach. 
It should be mentioned that a similar problem afflicts the instanton liquid model as well (when 
choosing size and density as quoted above). 

Our final goal is to complement the analytical treatment of the dyon model presented in ref. [l] 
by numerical simulations of the proposed moduli space metric together with interactions stem- 
ming from the gluonic fluctuation determinant, the Faddeev-Popov determinant and the action 
itself. The preparatory studies to be presented here will demonstrate that the approximate 
moduli space metric from [l] violates the fundamental requirement of being positive deflnite in 
an overwhelming part of dyon configuration space, when explored at densities needed to match 
phenomenology. We conclude that this moduli space metric can only be used in a dyon model 
with other terms added, or if correlations between the dyons are built in that lead to the posi- 
tivity of the metric. A suitable simulation algorithm should guarantee that. This would result 
in complicated multi-dyon correlations and create a non-trivially related behavior of various 
observables and correlators. 

This paper is organized as follows. In section[2]we briefly recall the ingredients of the dyon model 
in the version of Diakonov and Petrov applied to SU(2) gauge theory, which is the starting point 
for our work. Section [3] deals with the spectral properties of the proposed moduli space metric. 
In section |4] we demonstrate that even a random model of dyons, i.e. a model without moduli 
space metric, induces confinement. Finally we conclude and give a brief outlook. 

2 Dyon ensembles a la Diakonov and Petrov 

2.1 Holonomy as external parameter 

We consider pure SU(2) gauge theory at finite temperature T, which, as usual, is implemented 
with periodic boundary conditions in the imaginary time direction with period [5 = 1/T. There 
are four kinds of dyons, two selfdual and two antiselfdual. The untraced Polyakov loop at spatial 
infinity, also called holonomy, is an element of the gauge group. It can be diagonalized every- 
where! to = exp(27rzL(j(T3) with eigenvalues e"''^'^*'^ and e"^'^*'^, where lo S [0, 1/2]. Similarly to 
the non-Abelian adjoint Higgs system this gives rise to complementary dyons that have actions 
proportional to 2uj and 1 — 2a; and that become static when well-separated. They are both 
selfdual or antiselfdual depending on the sign of their topological charge. 

Following Diakonov and Petrov [1], throughout this paper we focus on the case of maximally 

^To be more precise, it can be diagonalized everywhere except for the loci of the Dirac strings, if the system 
is not neutral. 
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non-trivial holonomy, oo = 1/4. This implies that all dyons have the same action. The holonomy 
is then traceless and matches the confinement condition {P) = 0, where P = Tr(r2)/2. We adopt 
the same restriction to purely selfdual systems, considering K dyons of the first kind and K 
dyons of the second kind, i.e. the total number of dyons \s no = 2K. These dyons carry equal 
electric and magnetic charges which can take the values ±1 . 



2.2 The approximate moduli space metric of multi-dyon configurations 

In [T] it has been attempted to construct an approximate multi-dyon moduli space metric valid 
for dyon separations d ^ l/vrT. The starting point for this construction were the analytically 
known moduli space metric of a single caloron [31], which is a pair of different kind dyons, and 
a corresponding approximation for pairs of same kind dyons at large distances. The integration 
over collective coordinates, which are the dyon positions x™" [f| with i = 1, . . . denoting the 
dyon index and m = 1,2 denoting first and second kind respectively, is then performed with the 
measure 

K 2 X 

\{J{d^xT)V^)- (1) 

i=l m=l ^ 

The approximate moduli space metric g is related to a matrix G, 

/ ^^1 1 \ 5™'"(1 — 5i •) 

V fc=l \k k=l,kjti "i.fc ^ "^.i %j 

(see appendix |A]) such that the determinants are related by 



Vdet{g) = det(G). (3) 

The notation c?™;" = |x™ — x"| is used for dyon separations, and dyon type indices m are 
considered as cyclic, i.e. for SU(2) holds (m = 3) = (m = 1). We omit the temperature T 
in most equations, which means that distances are measured in 1/T units. The temperature 
dependence can easily be restored by making this explicit. 



2.3 Parameters and their physical values 

Within the model proposed in [1] the relation between the string tension a and the critical 
temperature Tc of the confinement deconfinement phase transition can be derived analytically. 
It reads for SU(2) gauge theory: 

a = (SpTj , = (iSpT/TT^j , (4) 

or equivalently 



T. / 6 ^ 



0.883. (5) 



■^The phases of the dyons are irrelevant in this context. 
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Eqns. dH) are given in parametric form, where p is the three-dimensional density of dyons of 
each kind, i.e. p = K/V^, = riD/'^y'i with V3 denoting the spatial volume. Hence, pT can be 
interpreted as the four-dimensional density of dyons which should actually be viewed as the 
fundamental parameter of the model. 

In order to fix physical units, we set the scale by choosing (3 = l/T = l.OOfm = 1/198 MeV 
throughout the paper. This means, we are considering a certain temperature in the confining 
phase of SU(2) Yang- Mills theory. Setting the string tension to its "physical value" , taken here 
as cTphysicai = (440 MeV)^ = 4.99/fm^, eqns. <^ yield a three-dimensional density po = S.lO/fm"^ 
and Tc = 389 MeV. The dimensionless ratio Td^fo = 0.883 differs from that quoted in [T] by 
the factor 2^/^, which originates from the replacement o — > \/2a made in [T] in order to take 
the contribution of antiselfdual dyons into account. In general, we define d= (1/p)^/^, which 
is proportional to and of the same order of magnitude as the average nearest neighbor dyon 
separation in three-dimensional space at density p. In particular, Jq = (1/po)^''^ = 0.686 fm is, 
according to ([4]), the dyon separation reproducing the physical value of the string tension for 
our choice of temperature, T = 198 MeV. 

3 Spectral properties of the approximate multi-dyon moduli 
space metric 

We consider selfdual multi-dyon configurations, which are solutions of the classical Yang-Mills 
equations of motion and, therefore, local minima of the action functional of Yang-Mills the- 
ory. The corresponding surface of minimal action can be parameterized by a set of collective 
coordinates, the dyon positions x™ (which have already been introduced in section 12. 2p and 
phases which are unimportant for the following. This surface is embedded in a flat Euclidean 
space of infinite dimension spanned by the gauge field degrees of freedom AJ^(x). Consequently, 
the induced moduli space metric associated with these collective coordinates must he positive 
definite. 

The approximate multi-dyon moduli space metric g proposed in [T] and the corresponding de- 
terminant have been given in eqns. ([T|) to ([2]). For the investigation of positive definiteness 
one has to consider the eigenvalues of g. In appendix [A] we show that the number of negative 
eigenvalues of g is equal to four times the number of negative eigenvalues of G. Thus it suffices 
to study G: if G is not positive definite, the same holds for g, which means that it fails to satisfy 
the fundamental property of positive definiteness inherent to any moduli space metric. In such 
cases the weight factor det(G) associated with the corresponding dyon configuration cannot be 
interpreted in a physically meaningful way. 

Since G is only an approximation of a moduli space metric, the existence of loci of non-positive 
definiteness are not excluded. Still, one could hope that such cases of failure of the used approxi- 
mation are restricted to a small part of dyon configuration space, which might even vanish in the 
thermodynamic limit. Our numerical investigations, however, strongly indicate the opposite, i.e. 
that even at small dyon densities the percentage of configurations with positive definite G tends 
to 0, when the size of the system is increased. In the following we present a detailed study of 
the spectrum of G, particularly its positive definiteness, for various dyon numbers and densities. 
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3.1 Analytical considerations 



We start by considering the simple case of a pair of different kind dyons, i.e. a caloron. The 
moduli space metric is exactly known [31]. At maximally non-trivial holonomy it is given by 

- \ -2/d 2n + 2/d ) ' 
where d is the dyon separation. It is positive definite with eigenvalues 

Xi = 2tt , \2 = 2TT + A/d. (7) 
Hence, the weight factor det[Gd) is positive for arbitrary dyon positions. 

Considering, on the other hand, a pair of same kind dyons (using the approximate moduli space 
metric proposed in [1]) the signs in front of the 2/d terms are reversed, 

_ ( 27, -2/d +2/d \ 
~ V +2/^^ 2n-2/d J ' 

yielding eigenvalues 

Ai = 27r , A2 = 2-71 -4/d. (9) 

Here the weight factor det(Gs) is positive only for dyon separations d > 2/tt = 2/ttT. For 
our choice of temperature, T = 198 MeV, this corresponds to d > 0.635 fm, which is of the 
same order of magnitude as the average distance of same kind dyons dQ = 0.686 fm needed to 
reproduce the physical value of the string tension CTphysicai = 

(440MeV)2 (cf. section E^]). It has 
been argued that at the critical d the distance ceases to be a good collective coordinate (see [l] 
and references therein) as the dyons strongly overlap. 

Of course, such a pair of same kind dyons has to be complemented by another pair of the other 
kind to arrive at the neutral situation that we are going to investigate in general (in other words 
two-caloron solutions need to be studied, of which a few are known [321 133j). Then the moduli 
space metric becomes a 4 x 4 matrix and the conditions, under which some of its eigenvalues are 
negative, are not that obvious. In section [3^2] and [3.31 we perform a detailed numerical analysis 
of the dependence of the spectrum of G on the dyon number and density both for randomly 
sampled dyon positions and for dyon positions distributed according to | det(G)|. The examples 
presented here already indicate that the approximate moduli space metric G might fail to fulfill 
the fundamental requirement of positive definiteness. 

In the following we show that for fixed dyon number n£) the percentage of configurations with 
non-positive definite G becomes larger, when the dyon density is increased. The matrix G can 
be written in the form 

G = 2Tlln^^nr, + D{^J , (10) 

where the elements of D are linear combinations of inverse dyon distances 1/d™^" (cf. eqn. 
D is, therefore, inversely homogeneous in the dyon positions, i.e. 

D{^Ja) = aDi^l,). (11) 
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We start by considering an arbitrary dyon configuration {x^} at density p fulfilling 

Tr(G) < 27rnz). (12) 

For large dyon numbers n^), roughly half of the configurations will satisfy this requirement. This 
is so, because Tv{D) is a sum of {ud/'^Y positive and nD/2 x (n£)/2 — 1) negative terms, i.e. 
it fluctuates approximately around zero, and because a traceless D gives exactly the bound of 
eqn. (jl2p . Obviously, such a matrix G has at least one eigenvalue smaller than 27r, but might of 
course be positive definite. 

Now we scale all dyon locations x™ by 1/a < 1, i.e. {x™} {x'™} = {x™/a}, which amounts 
to increasing the dyon density by the factor a'^, i.e. p ^ p' = a^p. Because of the inverse 
homogeneity of D (eqn. (fTTj) ). the eigenvalues of the squeezed configuration, A^, are related to 
the eigenvalues of the original configuration, Aj, by multiplying their difference to 2tt by a, i.e. 

\'j = 2tt + a{Xj - 2tt) . (13) 

Consequently, for any eigenvalue Xj < 2tt one can choose a such that A^- < 0. In other words, 
by rescaling any such generic configuration will give rise to a configuration for which G is not 
positive definite. 

This argument can be transfered to the average spectral density of G. The latter is obtained 
by averaging over randomly and uniformly chosen dyon positions inside a cubic spatial volume. 
Comparing spectral densities for dyon densities p and p' (at fixed dyon number no) one finds 
that the latter is stretched by the factor a around a "fixed center" at 2tt (cf. eqn. ()13p and 
also Figured] for numerical evidence, where the center at 27r is indicated by dashed lines). This 
illustrates that all eigenvalues smaller than 27r will eventually become negative, when the dyon 
density is increased. 



3.2 Multi-dyon configurations at fixed density: dependence of the spectrum 
of G on the dyon number 

In the following we consider dyon configurations with randomly and uniformly chosen positions 
inside a cubic spatial volume V3 = L^. At fixed density p we gradually increase the dyon 
number rir, (and consequently V3), to investigate the effect on various quantities characterizing 
the positive definiteness of G: 

• the average percentage R of negative eigenvalues of G, i? = {n^)/n£), where n_ denotes 
the number of negative eigenvalues for a given dyon configuration; 

• the probability distribution P~{n-) of the number of negative eigenvalues of G; 

• the spectral density of G (obtained by averaging over sampled dyon configurations). 

Results for two significantly different choices of the dyon density, p = 1/fm'^ and p = 1/125 fm^, 
are collected in Table [T] and Figure [TJ 

At first we discuss p = 1/fm^, a density smaller by more than a factor of 3 than the density 
Pq = 3.10/fm^, which is according to [1] needed to produce a quantitatively correct string 
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no = 100 




0.4 n^/uD 



0.4 n^/uD 




0.4 n^/riD 




0.04 



p = 1/125 W 
no = 200 



0.04 n^/no 



p = 1/125 W 
no = 400 



0.4 n^/riD 




0.04 n^/no 



1/125 W 
no = 800 



0.02 



0.04 



11- /no 



Figure 1: Histograms obtained from 100,000 independently chosen dyon configurations repre- 
senting the probabihty distribution P~{n-) of the number of negative eigenvalues (as function 
of n_/n£)). Left column: p = 1/fm'^. Right column: p = 1/125 fm'^. 
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p = 1/125 W 


no 


L in fm 


R in % 


p-{0) in % 


L in fm 


R in % 


p-(0) in % 


100 


3.68 


21.81(6) 


0.000(0) 


18.42 


0.457(5) 


64.53(35) 


200 


4.64 


25.04(3) 


0.000(0) 


23.21 


0.546(2) 


36.72(2) 


400 


5.85 


28.16(2) 


0.000(0) 


29.24 


0.684(1) 


10.22(3) 


800 


7.37 


30.94(1) 


0.000(0) 


36.84 


0.919(4) 


0.60(40) 



Table 1: The average percentage R of negative eigenvalues and the percentage of dyon config- 
urations with positive definite G, P~{n- = 0), for various dyon numbers no and two selected 
densities p (the averages have been computed from 100, 000 independently chosen configura- 
tions). 



tension, etc. (cf. also eqn. ([ID). The left columns of Table [T] and Figure [T] clearly show that for 
dyon numbers in the range of 100 < n£, < 800 the probability -P~(0) to find a configurations 
with positive definite G, is essentially zero. In other words, among the 100, 000 independently 
chosen dyon configurations there is not a single configuration, where the weight factor det(G) 
can be interpreted in a physically meaningful way, i.e. coming from a positive definite moduli 
space metric. Even worse, the average percentage of negative eigenvalues R increases, when we 
consider a larger number of dyons. This implies that a model based on the approximate moduli 
space metric G and enforcing its positive definiteness contains highly non-trivial constraints 
determining a small admissible subspace of configuration space, and that the restrictiveness of 
these constraints even increases in the thermodynamic limit. 

Since the construction of G in [1] is based on the assumption that all dyons are well separated, 
one could hope that reducing their density might cure the problem. To check this, we have 
repeated the analysis for a rather dilute ensemble with p = 1/125 fm^. Note that according to 
the model proposed in [1] this dyon density would yield a rather low value for the string tension, 
a w 0.05 X (Tphysicai (cf. eqn. Indeed, the average percentage of negative eigenvalues R is 

now significantly smaller, as can be seen in the right columns of Tabled] and Figure [TJ However, 
as before, this percentage increases, when we consider a larger number of dyons (keeping the 
same density), e.g. for n£) = 800 dyons less than 1% of all configurations have an associated 
matrix G, which is positive definite. This gives additional evidence that even for rather low 
densities the positive definiteness of the moduli space metric remains a very selective constraint. 

One could still hope that results obtained by an unrestricted integration over collective coor- 
dinates might evade this problem in the sense that the integrated weight associated with dyon 
configurations with positive definite G, W(j(n_ = 0), might be significantly larger than the in- 
tegrated weights corresponding to "unphysical sectors" , where G is not positive definite, namely 
Tyj(n_ > 1): 



/ (n n ^'^r) A„_(xr) I det(G(xr)) 

^ i=l m,=l ^ 

nn^'^r det(G(xn) 



Wain.) 



Z 



K 2 



=1 m=l 



(14) 
(15) 
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where 



1 if G(x™) has n_ negative eigenvalues 
otherwise 



(16) 



We investigate this possibihty by computing the integrated weights with the absolute value of 
det(G) and sorting them with respect to the number of negative eigenvalues n_. If Wq{0) 
happens to be much larger than VF(J(n_ > 1), then the positivity problem might be less severe 
in actual simulations. In ref. [1] this was expected from the strong repulsion of dyons at a zero 
of the metric determinant. 

Since det(G) exhibits strong fluctuations covering many orders of magnitude, we have evaluated 
via Metropolis sampling writing |det(G)| = exp(ln | det(G)|). Results at dyon density 
p = 1/fm^ are shown in Figure EJ It is clearly visible, that the integrated weight Wq{0) of the 
physically meaningful sector without negative eigenvalues is negligible compared to the absolute 
weight of all unphysical sectors, X]„_>i VF^ (n_). Moreover, one can see that for dyon numbers 
in the range 50 < n£) < 200 physical observables are dominated by dyon configurations, where G 
has around 80% to 90% negative eigenvalues, and where follows a smooth bell-shaped 

curve. 
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no = 100 
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20r 



15 



10 



no = 200 
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Figure 2: Integrated weights Wq {n^ 
n£) and a fixed density p = 1/fm^. 



(shown as functions of n-/n£)) for various dyon numbers 



Taking the full determinant into account, i.e. taking det(G) instead of |det(G)|, one has to 
multiply iy^(n_) by an alternating sign (+ if n_ is even, — if n_ is odd). Then roughly 
half of the dyon configurations have associated negative weights, when one computes ensemble 
averages ignoring the requirement of positive definiteness. Similar to the case with randomly 
chosen dyon positions, the average percentage R of negative eigenvalues of G increases, when 
the dyon number is increased at fixed density. 

We conclude that results based on the approximate moduli space metric G by performing an 
unrestricted integration over collective coordinates (without control over positive definiteness) 
are physically meaningless. 
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3.3 Multi-dyon configurations at fixed dyon number: dependence of the 
spectrum of G on the density 

To investigate the dependence of the spectrum of G on the dyon density p, we proceed in the 
same way as in section 13. 2| this time keeping the dyon number n£, fixed, while varying the 
spatial volume V3. 

For randomly and uniformly chosen dyon positions, our results show that the average percentage 
of negative eigenvalues R becomes larger, when the density is increased (cf. Tableland Figure^]). 
This is in agreement with the analytical argument given in section [3Tl This result was expected, 
since the moduli space metric G is an approximation valid only for large dyon separations. 
Even when considering rather dilute ensembles (p = 1/125 fm'^), still less then 40% of all dyon 
configurations for n£, = 200 and less than 1% of all dyon configurations for riu = 800 possess a 
positive definite G. 





no = 200 


UD = 800 


p 


L in fm 


R in % 


p-(0) in % 


L in fm 


R in % 


p-(0) in % 


1/125 fm^ 


23.21 


0.546(2) 


36.72(2) 


36.84 


0.919(4) 


0.60(40) 


1/64 fm^ 


18.57 


1.225(8) 


13.28(12) 


29.47 


2.305(8) 


0.006(4) 


1/27 fm^ 


13.92 


3.241(3) 


1.183(12) 


22.10 


5.997(7) 


0.000(0) 


1/8 fm-'^ 


9.28 


9.382(18) 


0.003(1) 


14.74 


14.60(2) 


0.000(0) 


l/fm^ 


4.64 


25.04(3) 


0.000(0) 


7.37 


30.94(1) 


0.000(0) 


po = S.lO/fm^ 


3.18 


33.08(2) 


0.000(0) 


5.05 


38.10(1) 


0.000(0) 



Table 2: The average percentage R of negative eigenvalues and the percentage of dyon configu- 
rations with positive definite G, P'~{n- = 0), for two dyon numbers no and various values of the 
density p (the averages have been computed from 100, 000 independently chosen configurations). 



To complement the picture, we show in Figured the spectral density of G both for n£, = 200 
and for no = 800. From these plots one can clearly see that large dyon numbers or densities 
inevitably yield configurations, where the associated matrices G are not positive definite. The 
scaling of the spectrum predicted in section 13.11 is nicely confirmed. 

We have also studied the dependence of the integrated weights W(^(n_) (cf. eqn. (fH|) ) on the 
dyon density p at fixed dyon number n£) = 200. Results are shown in Figure [5j Again one 
can see that for p > 1/8 fm^ observables are dominated by dyon configurations with associated 
matrices G, which have between 80% and 90% negative eigenvalues. Moreover, the maximum of 
the integrated weights is shifted towards larger numbers of negative eigenvalues, when the dyon 
density is increased towards po that is considered a realistic value. 
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Figure 3: Histograms obtained from 100,000 independently chosen dyon configurations repre- 
senting the probabihty distribution P~{n-) of the number of negative eigenvalues (as function 
of n^/no). Left column: n£, = 200. Right column: no = 800. 
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Figure 4: Histograms obtained from 100,000 independently chosen dyon configurations repre- 
senting the spectral density of the approximate moduli space metric G. Left column: nu = 200. 
Right column: ud = 800. 
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Figure 5: Integrated weights T4^(^(n„), shown as functions of n-/n£,, as obtained for a number 
n£) = 200 of dyons and various densities. 

4 Confinement from a random dyon gas 



In the fohowing we consider dyon ensembles without moduh space metric (or other interac- 
tions), i.e. we perform a uniform samphng of dyon positions. Although this is a rather extreme 
simplification of a dyon model (see also [23]), such ensembles might be worth studying, be- 
cause confinement still persists, as we will shortly demonstrate. This investigation might also 
be helpful to better understand the impact of a moduli space metric on ensembles of dyons. 

We consider ud = 1,600 dyons at the same densities as in section [3l We compute the Polyakov 
loop correlator yielding the free energy of a pair of static charges: 

Fqq{R) = -rin(p(x)pt(y)^ , i? = |x-y|. (17) 
The Polyakov loops 

P(x) = ^Tr(^exp (^f^ ^ dxo^o(x))) = cos (4(x)/2t) . (18) 

can be evaluated analytically, because the gauge field is treated in the Abelian far field limit 
(see appendix [B]) . In the algebraic gauge it is given by 

ao(x;g) = - , ai(x;g) = , . , al{:>c;q) = + /^^ , (19) 

r r(r — z) r{r — z) 

where x = (x,y,z) and r = |x|. All other gauge field components are zero. The multi-dyon 
configurations we use are linear superpositions of gauge potentials (I19p determined by dyon 
positions x^ and charges q"^, 

K 2 

A;1(x) = <^«35^ovrT + J] J]a«(x-xr;g™), (20) 

i=l m,=l 
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where g'" = ±1 for the m = 1,2 dyons plus an additional term in generating (in periodic 
gauge) the non-trivial holonomy 0, = exp(z7r/2cj3) sufficiently far away from all dyons. We 
regularize the singularities at the dyon centers by factors 1 — exp(— r/e) with e = 0.25 fm0. 

As is shown in Figure E] the free energy rises linearly up to w 1.5 fm for p = 3.10/fm'^ and 
up to ~ 2.7 fm for p = 1/fm^, where statistical noise starts to dominate and up to ~ 4fm for 
P < 1/8 fm^, which is the maximum separation considered. 




0.5 1 1.5 2 2.5 1 2 3 4 

separation R in fm separation R in fm 



Figure 6: The free energy of a pair of static charges Fqq as a function of the separation R for 
various dyon densities p. 

The corresponding string tensions are summarized in Table [3l We conclude that even non- 
interacting dyons give rise to confinement, when considering the Polyakov loop correlator. Quan- 
titatively there is a strong dependence of the string tension a on the dyon density p. Note that 
for dyon density p = po = 3.10/fm^ (which is according to the model proposed in [l] that density 
yielding approximate agreement with lattice results) the extracted value of the string tension 
even slightly overshoots its physical value: a/a^^^^^'^^^ ~ 1.4. 



p 


L in fm 


fj/o-Physical 


1/125 fm^ 


46.42 


0.015(1) 


1/64 fm^ 


37.13 


0.020(1) 


1/27 fm^ 


27.85 


0.037(1) 


1/8 fm^ 


18.57 


0.074(1) 


l/fm^ 


9.28 


0.466(13) 


Po = 3.10/fm^ 


6.37 


1.434(19) 



Table 3: String tensions in units of the physical string tension cjphysical = (440 MeV)^ for various 
dyon densities p. 

On the other hand it is hardly surprising that confinement is already present on such a simple 
level. Due to the long range nature of the dyon far fields, the models studied in [T] and also in this 
paper exhibit certain similarities to ensembles of regular gauge instantons and merons [34^ [35l I36j 

*We have checked that numerical results remain essentially unaltered, when e is further decreased. 
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and to the pseudoparticle approach [37 1 138 1 [39 11401141] . foi^ which it is well known that confinement 
and the long range nature of the "gauge field building blocks" are closely related (in particular 
cf . ! where it has been demonstrated that even a random ensemble of merons yields a linearly 
rising static potential, as well as [381139). where the relation between the long range nature of the 
gauge field building blocks and confinement has been established). In a semiclassical approach 
at finite temperature, however, dyons are rather the natural building blocks. 

5 Conclusions 

The results obtained in |lj are based on an approximation of the multi-dyon moduli space metric, 
as shown in eqn. ([2]). This metric G - as we have demonstrated above - is positive definite only 
at small dyon numbers nD and small dyon densities p. However, it is used at values of n^j and 
yO, where it does not satisfy the fundamental requirement of positive definiteness. 

In detail our findings are the following: 

• At dyon numbers 100 < tid < 800 and typical densities {p ~ 1/fm'^) practically all dyon 
configurations correspond to matrices G, which are not positive definite. 

• Roughly half of the dyon configurations have odd numbers of negative eigenvalues. This 
implies that every second dyon configuration receives a negative weight factor, when the 
metric determinant is taken into account as weight of an unrestricted integration over 
collective coordinates. 

• All attempts to approach the thermodynamic limit by increasing the dyon number, while 
keeping their density fixed, have lead to more severe violations of positive definiteness. 

• Decreasing the temperature (throughout this work we have used (3 = 1/T = 1.00 fm), 
while keeping the physical values of the string tension and the critical temperature fixed 
(according to eqn. ([4])), amounts to increasing the three-dimensional density. As we have 
demonstrated this makes the situation even worse. 

• Dyons generically induce confinement; already dyon ensembles with randomly chosen po- 
sitions do so. 

We expect that the problems encountered for SU(2) are generic also for higher gauge groups. 

We consider these findings a challenge. Since positive definiteness is a fundamental property 
of any moduli space metric, it seems doubtful that results without restriction to the "positive 
definite subset" of dyon configurations can be interpreted in a physically meaningful way. If 
the metric determinant is taken into account as weight, averages of physical observables are 
computed from alternating sums over the number of negative eigenvalues of G. We conclude 
that, in order to obtain physically meaningful results, it is necessary to either modify the dyon 
model in a way that the (corrected) weight factor is always positive, or to restrict the integration 
over dyon positions to those parts of configuration space, where G is positive definite. 

In a subsequent publication we plan to present numerical simulations of dyon models with 
interactions derived from the approximate moduli space metric G, which do not suffer from 
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"unphysical sign problems" . One appealing possibility is the following integration over collective 
coordinates: 

^i=lm=l ^ 

Since the constraints imposed on the dyon coordinates by this measure are highly non- 
trivial, it seems unlikely that this model can be treated analytically. Therefore, we are currently 
developing efficient Monte-Carlo algorithms, which are respecting the positive definiteness of 
G, i.e. which are designed to avoid "forbidden" multi-dyon configurations. In this respect our 
finding that already dyons with randomly chosen positions generate confinement is encouraging 
in the sense that models based on dyons seem to capture the relevant degrees of freedom of 
SU(2) Yang-Mills theory. 



A Relation between the spectrum of G and the spectrum of g 

In the following we show that the number n-{g) of negative eigenvalues of g is four times the 
number (G) of negative eigenvalues of G. 

The relation between G and g is 



/ G + W^G-^W^ 
W^G-'W, 
WjG-^W^ 

/ 1 I^J \ 
1 
1 

V 1 / 



W^G"^Wy 
G + W^G~^Wy 

wTg-^w„ 



y 

wT 



G~^W, 



( G 



V 




G 





y 


G 




w]:g~^w, wJg-^ \ 



W^G-^W, 
G + W^G-^W^ 
G-^W, 



















1 





wTg-' 



G- 





1 



\ 




1 / 



(22) 



--W 



cf. eqn. (20) in [T]. 

Obviously, G has 4n_(G) negative eigenvalues. We denote the eigenvectors of G by v^'^'^^ 
and the eigenvalues by A^'^'^-*, i.e. Gv^'''^^ = A^'^V^'^'^^. -|- denotes a positive eigenvalue, i.e. 

> 0, and — denotes a negative eigenvalue, i.e. A^'^'"-' < 0. Because G is symmetric, the 
eigenvectors can be chosen orthonormal. Moreover, all eigenvalues are real and all eigenvectors 
can be chosen real. 

With arbitrary real coefficients a{k) one has 



> ^a(A;)2A('='-) = ( ^ a(A;)v('='")) g( Jj^ a(A;')v(^''-)) 



(23) 
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The 4n_(G) dimensional subspace 

S = W-^^a{k)^r^'''-^ (24) 
k 

can be expanded in terms of the orthonormal and real eigenvectors x^^'^^ of g (g is also sym- 
metric), i.e. 

S = ^/?(A:,-)x('='-)+^/3(A;,+)x('='+). (25) 

k k 

If g would have less than 4n_(G) negative eigenvalues, one could chose a combination of a{k) 
such that P{k,—) = 0. This, however, would be a contradiction to ([23]) . Therefore, g has at 
least 4n_(G) negative eigenvalues, i.e. n^{g) > 4n_(G). 

Analogously one can show that n^{g) > 4n+(G). This proves that n-{g) = 4n_(G). 



B The dyon gauge field in the far field limit 

The gauge field of a single dyon can be obtained by considering the gauge field of a caloron at 
maximal holonomy in the limit of infinite dyon separation. 

When the distance to the dyon center r = (x^ + y^ + z^)^/^ is large, the gauge field Af^ = A'^^a"' /2 
is Abelian: 

4 = 0,4 = 0,4 = 5^^TTT + qf^,d,Hct>), (26) 
where 

1 



r — z 



(27) 



and f)'^^ = eapiu — (5a^t<^0i/ + ^au^Oii is the t'Hooft symbol. The charge q is either +1 (dyons of the 
first kind) or —1 (dyons of the second kind). The coordinate system has been chosen such, that 
the singular Dirac string points in positive z-direction. 

The non-zero components of the gauge field read 

4 = vrT + ^ , 4 = ' ^2 = +T^- (28) 

r r[r — z) ryr — z) 

With the definitions Ej = Fqj, Bj = —ejkiFki/2 and F^j^^ = dfj^A^ — d,yA'^ the corresponding 
electric and magnetic monopole fields are given by 

E = 5 , B = 5. (29) 
Throughout this paper we approximate dyons always by these Abelian far fields. 
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